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tive quantum number for a relative velocity of 10 km s are also given together 
with velocity exponents A, if they exist, on the assumption that the cross section 
varies with velocity as v~ x . A discussion of our results is presented. 
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Abstract. The theory of collisional depolarization of spectral lines by atomic hy- 
drogen (Derouich et al. 120031 is extended to d (1=2) atomic levels. Depolarization 
rates, polarization and population transfer rates are calculated and results are 
given as a function of the temperature. Cross sections as a function of the effec- 



1. Introduction 

In stellar atmospheres, isotropic collisions with the particles of the medium are responsi- 
ble for a part of the broadening of spectral lines. If a polarizing effect creates an atomic 
polarization of the lines, they can decrease or even destroy this atomic polarization. 
Spectral lines observed at the solar limb are linearly polarized by anisotropic scatter- 
ing of the incident solar radiation. Isotropic collisions with the particles of the medium 
depolarize the lines. In the solar photosphere and the low chromosphere, depolarizing 
collisions are dominated by isotropic collisions with hydrogen atoms. When these depo- 
larizing collisions dominate over any other radiative (or collisional) polarizing effects, the 
atomic levels become depolarized. Therefore depolarization rates, polarization and 
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population transfer rates by collisions with hydrogen are needed in order to interpret 
the observed polarization. 

In Derouich et al. (2003), referred to hereafter as Paper I, a semi-classical theory 
for depolarization of neutral atomic lines by collisions with atomic hydrogen has been 
developed and applied to p (I = 1) atomic states. This theory is an extension to collisional 
depolarization of the theory developed in the 1990's by O'Mara and collaborators (Anstcc 
HMH Anstee & O'Mara HMTllTMSI Anstee, O'Mara & Ross HMTI Barklem lTMSl Barklem 
& 0'Mara lT997l Barklem, O'Mara & Ross ll99fr>|) for line broadening by collisions with 
atomic hydrogen. The present paper is an extension of this theory to d (I = 2) atomic 
levels. In fact, this paper presents the first calculations of the depolarization and the 
collisional transfer rates for c?-atomic states. 

A great advantage of the present method is that calculations are not specific for 
a given perturbed atom. The transition matrix T may be calculated using Coulomb 
wavefunctions for the valence electron of the perturbed atom and is simply dependent 
on the effective principal number n* and the orbital angular momentum quantum number 
I (1=2 for d-atomic states). Therefore we can calculate depolarization cross sections, 
polarization and population transfer cross sections for any level of any atom, 
allowing computation for complex atoms. This is very useful for interpreting the so- 
called "second solar spectrum" (Stenfo & Keller IT997[1 . where depolarization rates for 
many levels are needed (Manso Sainzo & Landi Degl'Innocenti 2002). 

An extension of our theory to higher ^-values, aimed at a more complete interpreta- 
tion of the second solar spectrum, will be the subject of further papers. 



2. Description of the problem and summary of Paper I 

In our collision problem, the perturbed atoms collide with a bath of perturbing hydrogen 
atoms. The hydrogen atom is assumed to remain in its ground state during the colli- 
sion. The internal states of the perturbed atom are described by the spherical tensor 
components nlJ p k of the density matrix. Due to the isotropy of the collisions, the depo- 
larization rates, polarization and population transfer rates are g-independent. 
The contribution of the isotropic collisions to the statistical equilibrium equations is: 

J nlj Jz 

t—J^Ull = -D k (nlJ,T) nU p L Q + Y, D k (nlJ - nlJ',T) nU p k Q (1) 
+ quenching term, 

where D k (nlJ, T) is the collisional depolarization rate for the statistical tensor of rank k. 
Each level of total angular momentum J relaxes with 2J + 1 independent depolarization 
rates. In particular D°(nlJ,T) is the destruction rate of population, D 1 (nlJ,T) is the 
destruction rate of orientation (circular atomic polarization) and D 2 (nlJ,T) is the de- 
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struction rate of alignment which is of interest in the understanding of the second solar 
spectrum. 

We assume that inelastic collisions with hydrogen perturbers do not alter the total 
population of an atomic level (nlJ). This is the so-called no-quenching approximation. 
The no-quenching approximation implies that 

D°(nlJ, T)=^2 C(nlJ -> nlJ', T) = 0, (2) 

where C,(nlJ — ► nlJ',T)(J ^ J') is the inelastic collisional rate. The expression for the 
depolarization rate D k (nlJ, T) is given in Paper I. 

If the quenching must be taken into account, D k (nlJ — > nl,J',T) corresponds to 
collisional transfer of population (k — 0), orientation (k — 1) and alignment (k = 2). 
D k (nlJ — > nlJ',T) can be written as a linear combination of the collisional transition 
rates between the fine stucture sublevels ^(nlJMj — > nlJ'M'j,T)(J ^ J') (Sahal-Brchot 
IT$77|> : 



D k {nlJ -^nlJ',T) = (2k + l) ^ ^y+J'-Mj-M'j | | 

Mj.M'j 




((nlJMj -> nlJ'M'j,T). 



J' k J' 

-M'j M'j 

The expressions between parentheses denote 3j-coefficients (Messiah[1961 ). 
In particular one obtains 



D°(nlJ -> nlJ', T) = ^j ^L±±((nlJ -> nlJ' , T), (4) 
where 

C(nlJ^nlJ',T) = ^— ((nlJMj^nlJ'M'j^T), (5) 

Mj,M'j 

is the collisional transition rate between fine structure levels. £(aJMj — » aJM'j, T) is 
the collisional transition rate between the sublevels \aJMj) — » \aJMj). It can be written 
as a function of the local temperature T and the hydrogen perturber local density tir 
(Paper I): 

y* 00 y* 00 

((aJMj -> aJM'j, T) =n H / 2irbdb v f(v) dv \{aJMj\I - S(b, v)\aJMj)\ 2 , (6) 

Jo Jo 

where / is the unit matrix and T — I — S is the so-called transition matrix depending 
on the impact-parameter b and relative velocity v. Then the collisional depolarization 
rates and the collisional transfer rates can be expressed in terms of the 5*-matrix elements 
for the collision which are functionally dependent on the interaction energy of hydrogen 
in its ground state with the perturbed atom. The essential difference between various 
theoretical computations of the depolarization and collisional transfer rates is in the 
method employed to determine the interaction energy and to determine the S-matrix. 
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3. Method 

Consider a collision between a perturbed atom in a I = 2 state and hydrogen in its ground 
state Is. The perturbed atom is described by an optical electron outside a positively 
charged core. The effective principal number is given by n* = [2(_Eoo — Eni)]^ 1 / 2 . E n i 
is the energy of the state of the valence electron and Eoo is the appropriate series limit 
for the parent configuration of the perturbed atom state. As in the p atomic levels 
case (Paper I), depolarizing collisions are due essentially to intermediate-range isotropic 
interactions between radiating and hydrogen atoms. The so-called Rayleigh-Schrdinger- 
Unsld (RSU) interaction potential, used in this work, is of semi-classical nature and is 
totally electrostatic. For more details we refer to Paper I and to the ABO papers (Anstee 
HM2l Anstee & O'Mara HMTllTMSI Anstee, O'Mara & Ross lH)9Tl BarklemQUO Barklem 
& O'Mara lTMTI Barklem, O'Mara & Ross HTM)) . We use the notation of Paper I in the 
following analysis. 

Inelastic collisions which leave the radiating atom in a final state (n'l 1 ) different from 
the initial one (nl) are neglected. Thus we consider only the (21 + 1) states of the 
subspace (nl). Since potentials are computed in the rotating frame (Paper I), which 
is obtained from the fixed laboratory frame by means of the geometrical rotation R 

fj §)j the interaction potential is diagonal. The (21 + 1) RSU potential elements 
calculated here are (Anstee & 0'Mara|l992) 

V e ff, mi = ( H ){100\(A){nlmi\V e ff\nlmi) (A) |100) (H) (7) 

1 1 f +co 

= (Mi\V\Mi) - —((Mi\v\Mi)') 2 + — / P^ l I llmil (R,p 2 )dp 2 , 
t^p ^p Jo 

in atomic units which are used hereafter. Here |M;) = |n^m;)( J 4)|100)(^f), Ii\ m , \ are lengthy 

complicated analytic functions and E p = —4/9 is adopted. P n *i is the Coulomb radial 

wave function for the valence electron of the perturbed atom with quantum defect 5 = 

n — n* . 

The total wavefunction of the system (atom+perturber) is taken as the product 
of the wave function IV 7 ) (A) of the perturbed atom and that of hydrogen in its ground 
state |100) (ff) : 

|^) = ^)(A)|100) (H) , (8) 

and is developed over the basis formed by the products of the separated atoms states 
|Afj>: 

|^))=$> M( (*)e- iB V|M;) (9 ) 

Mi 

where E Q Mi is the eigenenergy of the system made-up of the two isolated atoms. 

The semi-classical coupled linear differential equations are obtained from the time- 
dependent Schrodinger equation. The transformation from the rotating to the fixed lab- 
oratory frame is included following Roueff l|1974J) . For d-states, the coupled differential 
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equations become explicitly (Barklem & O'Mara 1997): 

.ao_2(t) = l _ a _ 2mVefffl + Weffi + Veff2) + J_ ao{t)e -2i0 { _ Vefffi + Veff2) 

+ l -a 2 {t)e^{W efffl - iV effA + V ef/t2 ) 
i ^T^ = \a-i{t)(V efL1 +V efL2 ) + \ ai {t)e-^{-V efL1 +V efL2 ) 
i^ = ^=«- 2 We^(-%, + V eff<2 ) + ±=a 2 {t)e-^{-V eff>0 + V eff<2 ) (10) 
+ \ao(t)(V eff>0 + 3V eff , 2 ) 

= \a^(t)e^(3V efffi - W effA + V eff , 2 ) + ^=a (t)e^(-V eff>0 + V eff<2 ) 

+ -^a 2 (t)(3V efffi + iV efL1 + V efL2 ) 

Having the interaction potential V e ff, after integration of these equations, we obtain the 
transition matrix elements in the \nlmi) basis for a given velocity and impact parameter. 
The T-matrix elements in the \nlJMj) basis, which are needed for the depolarization 
and collisional transfer rates calculations, are obtained from equation (21) of Paper I. 

In the irreducible tensorial operator basis, the angular average over all possible 
directions of the collision plane of the depolarization transition probability is given in 
Paper I. That of the collisional transfer transition probability is given by: 

(P k (nlJ -> nlf \b,v))av = (nU fJ-\T\nlJ' fi')(nlJ u\T\nlJ' u')* (11) 

X 

Owing to the selection rules for the 3j-coefncients, the summation over x is reduced to 
a single term, since \ = — W ~ A*') = — ( v ~ I- 1 )- 

The depolarization rates D k (nlJ, T), and the polarization and population trans- 
fer rates D k {nlJ — > nU' \T) follow from integration over the impact parameters and 
the velocities with a Maxwellian distribution (for more details see Paper I). 




4. Results 

In most cases, the behaviour of the cross sections with the relative velocity v obeys a 
power law of the form: 

a k (nlJ -> nlJ',v)(J = J' and J ^ J') = a k (nlJ -> nW, V() ){—)~ xk{nU ^ nlJ '\ (12) 

where vq is a typical velocity where the cross section is calculated (10 km s" 1 ). In certain 
cases here, such behaviour was not obeyed (the cross section showed oscillations with 
velocity). 



6 M. Derouich et al: Collisional depolarization and transfer rates. 



n* 


a 2 (n2|) 


a 2 (n22) 


a 2 (n2|) 


cr°(n2f -> n2|) 


cr 2 (n2| -» ri2|) 


2.5 


283 


507 


342 


278 


60 


2.6 


313 


566 


380 


311 


65 


2.7 


351 


633 


424 


351 


73 


2.8 


394 


717 


477 


397 


77 


2.9 


443 


808 


537 


451 


87 


3 


496 


900 


598 


512 


104 


3.1 


553 


1008 


666 


578 


116 


3.2 


628 


1109 


727 


647 


139 


3.3 


711 


1210 


800 


711 


163 


3.4 


812 


1348 


897 


783 


179 


3.5 


905 


1470 


964 


847 


190 


3.6 


1028 


1621 


1039 


920 


196 


3.7 


1108 


1763 


1106 


961 


183 


3.8 


1202 


1910 


1185 


1004 


195 


3.9 


1299 


2066 


1260 


1079 


206 


4 


1294 


2191 


1367 


1117 


232 



Table 1. Variation of the cross sections, for the relative velocity of 10 km s 1 , with the 
effective principal number. Cross sections are in atomic units. 

Table ^ gives the various cross sections as function of n* for a relative velocity of 
10 km s _1 and the corresponding velocity exponents X k (nlJ — > nlJ'), if the exponential 
behaviour was obeyed, are tabulated in Table Then we can obtain the cross sections 
a k (nlJ — > nlJ' , v) for all velocities from Tables 1 and 2 using equation l|12f> . Tables^and 
12 can be interpolated for an appropriate n* associated to a given observed line in order 
to obtain the needed depolarization cross sections and collisional transfer cross- 
sections. After integration over velocities of these cross-sections, we obtain 
the depolarization rates and the collisional transfer rates of the line studied. 

For cross sections obeying equation 112(1 . as in Anstee & O'Mara l|19H2jl . the depo- 
larization and the collisional transfer rates can be expressed by: 

D k (nlJ -> nlJ', T){J = J' and J ^ J') = (l)(i> fc M^»^'))r(2 - \\ k {nl.J -» nlJ 1 )) 

7T 2 

x v a a k {nlJ -> nlJ 1 1 v ){-) 1 - xk{nlJ ^ nlJ '^ (13) 

We can generalize this relationship, which is specific to a particular atom owing to its 
mass, by assuming that /i = m# (fx and mjj are the reduced and hydrogen mass respec- 
tively) . This approximation introduces a very small error (Paper I) . 

As mentioned, in certain cases especially for transfer of linear polarization calcula- 
tions, the cross sections do not show exponential behaviour with velocity and so A is not 
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n* 


A 2 (n2§) 


A 2 (n22) 


A 2 (n2§) 


A°(n2§ -»n2§) 


A 2 (n2§ -> 7i2|) 


2.5 


0.207 


0.258 


0.268 


0.259 


0.172 


2.6 


0.223 


0.260 


0.269 


0.258 


0.184 


2.7 


0.222 


0.253 


0.268 


0.259 


0.233 


2.8 


0.200 


0.253 


0.271 


0.257 


0.209 


2.9 


0.214 


0.261 


0.270 


0.253 


0.177 


3 


0.210 


0.265 


0.280 


0.249 


0.148 


3.1 


0.185 


0.273 


0.282 


0.245 


- 


3.2 


0.168 


0.274 


0.284 


0.257 


- 


3.3 


0.153 


0.255 


0.280 


0.267 


0.194 


3.4 


0.137 


0.253 


0.285 


0.279 


0.238 


3.5 


0.140 


0.258 


0.304 


0.286 


0.306 


3.6 


0.131 


0.247 


0.315 


0.272 


0.394 


3.7 


0.175 


0.257 


0.330 


0.249 


0.446 


3.8 


0.179 


0.250 


0.303 


0.234 


0.479 


3.9 


0.208 


0.283 


0.322 


0.237 


0.399 


4 


0.172 


0.262 


0.296 


0.226 


0.448 



Table 2. Velocity exponents X k (nlJ — > nlJ')(J = J' and J ^ J') corresponding to the 
cross sections of Table H 

reported (Table|2|l. In these cases the linear polarization transfer rates must be computed 
by a numerical integration over computed cross sections obtained at different velocities. 

Figure n snows the alignment depolarization rates (k = 2) as a function of the local 
temperature T and n* for I = 2. The population transfer rates (k = 0) and the linear 
polarization transfer rates (k — 2) as a function of T and n* are displayed in Figs [2] and 
13 All these rates increase with the temperature. For a temperature T < 10000 K, the 
destruction rate of alignment D 2 {n 2 3/2)/n# < 6 x 10~ 14 rad. m 3 s -1 , D 2 (n 2 2)/n# < 
8 x 10~ 14 rad. m 3 s" 1 and D 2 (n 2 5/2)/n H < 5 x 10~ 14 rad. m 3 s -1 . The population 
transfer rate D°{n 2 3/2 -> n 2 5/2)/n H < 4 x 10 -14 rad. m 3 s" 1 and the linear 
polarization transfer rate D 2 (n 2 3/2 n 2 5/2)/n H < 8 x 10~ 15 rad. m 3 s _1 . These 
numerical values are given for n* < 4 which include most of the lines of interest for the 
second solar spectrum studies. The linear polarization transfer rates are smaller than the 
other rates. In fact equation J3J) shows that D 2 (n 23/2— >n25/2)isa linear combination 
of C(nlJMj — > nlJ'M'j,T). The coefficients of this linear combination have the sign of 
(3M 2 — J( J + 1 ) ) x (3M' 2 — J' ( J' + 1 ) ) . Therefore these coefficients are sometimes positive 
and sometimes negative. Consequently, due to the compensation between the different 
collisional rates C,(nlJMj — ■> nlJ'M'j,T), D 2 (n 2 3/2 n 2 5/2) is small compared with 
the other rates. For similar reasons circular polarization transfer rates (k = 1) are negative 
for ( J = 1/2 — * J' = 3/2) and (J = 3/2 — * J' = 5/2). This remark is in agreement with 
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the negative quantum chemistry orientation transfer rate obtained by Kerkeni ( 2002 ) for 
Na D lines (J = 1/2 — > J' = 3/2). We do not give our results concerning the circular 
polarization transfer because they are not needed for interpretation of the second solar 
spectrum . 






Fig. 1. Depolarization rates (k=2), per unit H-atom density, as functions of temperature 
T and n* . For I = 2, each figure: S = \ and J = §; S = and J = 2; S = \ and J = §. 
Depolarization rates are given in 10~ 14 rad. m 3 s -1 . 
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9 




Fig. 2. Population transfer rate (fc=0), per unit H-atom density, as a function of tem- 
perature T and n*.Z = 2,5=|,J=| and J' = |. Population transfer rate is given in 
1(T 14 rad. m 3 s" 1 . 




Fig. 3. Linear polarization transfer rate (k—2), per unit H-atom density, as a function 
of temperature T and n* . 1 — 2, S = |, J = | and J' = |. Linear polarization transfer 
rate is given in 10~ 15 rad. m 3 s _1 . 

5. Discussion 

The interactions of importance for the depolarization rates, and the polarization and 
population transfer rates calculations are the intermediate-range interactions. In the 
corresponding regions accurate interaction potentials are required. Examination of the 
RSU potential, used in this work, shows that long-range regions have the usual R~ 6 
behaviour but intermediate-regions have ~ R~ 10 behaviour. The Van der Waals potential 
is proportional to R~ 6 at all separations and so that gives a good description of the 
problem only at long-range separations. The Van der Waals potential underestimates 
the magnitude of the intermediate-range interactions. For this reason it can be seen why 
the calculations using the Van der Waals potential underestimate the depolarization and 
collisional transfer rates values. The results of Paper I showed that the RSU potential 
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gave depolarization rates which were in agreement (< 20 % for solar temperatures) with 
the quantum chemistry calculations. 

Unfortunately, there is neither experimental nor quantum chemistry depolarization 
and collisional transfer rates for d states (at least to the authors' knowledge) to compare 
with. The main differences between the RSU potentials and those from quantum chem- 
istry, which are more accurate, occur at short-range interactions. We verified that, as for 
the p states calculations in Paper I, these close collisions do not influence the computed 
depolarization and collisional transfer rates for the d states. We expect that a rather 
good agreement with a full quantum mechanical treatment (difference less than 20 %) 
would also occur for our present d states results. 

6. Conclusion 

The problem is to determine a great number of depolarization and collisional transfer 
rates by isotropic collisions with H-atoms. Our method was presented, tested and used 
with success for p states in Paper I. In the present paper we have given depolarization 
and collisional transfer rates corresponding to d states. These results are needed to model 
the formation of observed lines, and thus interpret the observations in terms of the solar 
magnetic field. The need is particularly strong for data for heavy complex atoms which 
are inaccessible to the quantum chemistry approach. An extension to / atomic states 
(I = 3) is a further interesting step in view of an extrapolation for / > 3 states. Such 
an extension to higher Z-values would be useful for a global interpretation of the "second 
solar spectrum" . This work is in progress. Adaptation and application of our theory to 
the determination of the depolarization and collisional transfer rates of singly ionized 
atoms lines by collisions with H-atoms will also be the subject of further papers. 
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